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We point out that in the early universe, for temperatures in the approximate interval 175-80 
MeV (after the quark-gluon plasma), pions carried a large share of the entropy and supported the 
largest inhomogeneities. Thus, we examine the production of entropy in a pion gas, particularizing 
to inhomogeneities of the temperature, for which we benefit from the known thermal conductivity. 
We finally put that entropy produced in relaxing such thermal inhomogeneities in the broad context 
of this relatively unexplored phase of early-universe cosmology. 


I. INTRODUCTION 

The hadron and lepton phase of early-universe cosmol¬ 
ogy, spanning a temperature range between about 175 
MeV and 1 MeV (between the quark-gluon plasma and 
cosmonucleosynthesis), has received only moderate at¬ 
tention in the literature [T], in spite of it being very rich 
in terms of the number of particles and interactions there 
present, and the underlying physics being relatively well 
known. This is probably because the only relic particles 
left from that era in the universe expansion form the cos¬ 
mic neutrino background [5] that there is at present no 
hope to detect. In spite of the dearth of direct messen¬ 
gers from that era, it is important to pursue its study for 
future precision work in cosmology. 

Particularly, there are few studies of the earlier part 
of the interval, just after exiting the quark-gluon plasma 
around 175-150 MeV [3], when a significant fraction of 
the universe’s entropy is carried by strongly interacting 
particles such as pions; only as the temperature drops 
below about 100 MeV, their decays ttq —>■ 77 , tt~ —>■ 
etc. add this entropy to that carried by leptons and 
photons. 

Much information on the pion phase is available from 
theoretical studies pertaining to the field of Relativistic 
Heavy Ion Collisions. Particularly, transport coefficients 
have been well calculated in recent years [IHT] and can 
be applied to early-universe physics. This is our focuse 
in the present work. The particular problem that we will 
address is entropy production. Though most treatments 
assume that the universe’s expansion is adiabatic and al¬ 
ways at equilibrium, this is just the simplest hypothesis 
and one may fancy consider separations from that equi¬ 
librium. 

One can argue that the rates of the particle-physics 
processes characterized in the Standard Model are larger 
than the expansion Hubble factor H = d/a, so the hy¬ 
pothesis of chemical and thermal equilibrium is reason¬ 
able, and the universe expands and cools down adiabat- 
ically. We of course concur with the analysis. But one 
cannot discard large past fluctuations in temperature or 
other quantities that have not survived to our days pre¬ 
cisely because of the large equilibration rates damping 


them. So there is always a level of hypothesis involved. 

What is solid information is that the fluctuations in the 
Cosmic Microwave Background (CMB) are measured and 
found small (VT/T < 10“®). So one can opt for evolving 
large initial-state inhomogeneities so they are this small 
at the time of recombination, or for considering inhomo¬ 
geneities that are so small in size as to evade observation 
in the CMB (cosmological versus microscopic inhomo¬ 
geneities). Further, since we will consider a radiation- 
dominated epoch, no structure-formation process is in¬ 
volved [8j. 

The largest contribution to the total entropy at ze¬ 
roth order during most of the time, is due to the 
(quasi)massless species (photons, neutrinos and elec¬ 
trons), as will be reminded below. Yet in a hot gas, since 
transport phenomena are diffusive, the typical transport 
coefficient (to which entropy production and the relax¬ 
ation rate will be proportional), drops with the inverse 
cross-section. The case in point for this study is the ther¬ 
mal conductivity, k (x 1/a. This means that the largest 
inhomogeneities at a given stage of the universe evolu¬ 
tion will be found in the gaseous subsystem which, being 
relativistic, is affected by the largest cross-sections. 

In the particle phase with components that are pho¬ 
tons, leptons, and pions, the largest entropy production is 
thus likely to take place in the pion gas. Heavier hadrons 
are barely present already shortly after the decay of the 
quark-gluon plasma, for example the kaon multiplicity [9] 
is down by at least an order of magnitude respect to the 
pion multiplicity. Therefore, though in principle kaon in¬ 
homogeneities can diffuse and produce entropy, we will 
ignore the phenomenon altogether. 

This is because their cross-sections are dictated by 
the strong QCD interactions and are in the 10-milibarn 
range, way larger than (Debye-screened, electromag¬ 
netic) lepton interactions. Even letting aside inhomo¬ 
geneities, for r > 80 MeV, pions actually carry a larger 
portion of the total homogeneous-gas entropy than pho¬ 
tons (though not larger than that of leptons) because 
of their multiplicity, as will be shown below in figure 
(bottom plot). Thus, there are two reasons to explore 
entropy and entropy production in the pion gas itself. 
In the high temperature end just after hadronization of 
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the quark-gluon plasma, pions are large carriers of en¬ 
tropy. And second, they are the ones that can support 
the largest inhomogeneities, if any are present, because 
they are the ones opposing diffusion most. 

With this motivation, our concrete study will be 
to address the relaxation of a thermal inhomogeneity 
at temperature T + ST towards the surrounding en¬ 
vironment value of T, ignoring other quantities that 
may separate from equilibrium such as momentum 
distributions or chemical inhomogeneities. Then we will 
calculate the subsequent entropy production to have 
reference values that may be useful in future studies. 


II. ENTROPY IN THE HOMOGENEOUS 
FRIEDMANN-ROBERTSON-WALKER 
COSMOLOGY 

A. System of equations for universe evolution 

In this section we quickly review the standard statisti¬ 
cal physics in the spatially flat (k = 0) [TO] homogeneous 
cosmos that serves as background for later study of inho¬ 
mogeneities. In this case the two independent Einstein 
equations give rise to the the Friedmann equation: 



for the evolution of the expansion parameter a(t) from 
the Friedmann equation and the balance equation 

^ = -^(p+P). ( 2 ) 

at a 

where P is the total pressure and the energy density p 
is the sum of the partial energy densities for the various 
species 

p = Pj + Pv,^ + Pe± + PiJ.,11 + ,7r0 + Pn,N + • • ■ (3) 

In table H] we summarize the main interaction channels 
in the temperature range we are discussing, 1 MeV < 
T <175 MeV. In particular, for entropy considerations, 
nucleons (already non-relativistic) and dark matter are 
not important. Pions and muons behave as radiation in 
the upper end of the temperature range. 

This can be seen for each species i, with degeneracy 
Pi, contributing 

= J > ( 4 ) 

because in thermal equilibrium, the function fi (usual 
Fermi-Dirac or Bose-Einstein distribution) 


suppresses the contribution by Here we have 

considered the more general case of local instantaneous 
thermodynamic equilibrium which will be useful later. 
As usual, p°‘ and t/“ are the components of the particle 
four-momentum p°‘ = (E, p) and the fluid four-velocity 
C7“ = 7 v(l, V) respectively, with yv = (1 — In 

a comoving frame, C/“ = (1,0,0, 0) and thus PaU°‘ = E. 
We have determined the chemical potentials pi from the 
current abundances and the scale factor a(t)/a(today). 
As they are tiny we do not quote them here. 
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Table I: Main processes in the temperature intervai 175 MeV— 
IMeV. Except when specified, the reactions can be written 
for ail charge combinations, e.g. tt denoting either of tt"*", 7 r“ 
or 7 r°, and p both muon and antimuon. (We omit additional 
reactions with much smaller branching fractions, such as the 
Dalitz decay tt® — >■ 77 * — >■ e^e “7 at 0(1%) [TT].) In this 
article we focus on the temperature range where tttt tttt 
dominates transport. 


Summing Eq. Q over species yields piT) from which 
the temperature evolution can be extracted as 


dT 

dt 


3a 

a 


ip + P) 


dT 

dp 


( 6 ) 


Eq. Q and ^ can be solved numerically by using for 
example the Runge-Kutta algorithm. The energy density 
is computed from the numerical integration of Eq. Q at 
each temperature, and the pressure is similarly obtained 
from the spatial trace of the energy-momentum tensor 
= 3P, which results in a sum over partial pres¬ 
sures of all species, 


P = ■ (7) 

i i 


Density and pressure decrease monotonically with t, 
while the scale factor a{t) increases monotonically; any 
of them may be used as a clock for further computations. 
We will set as origin of time the exit from the QGP at the 
top of the temperature interval, 0 = tT=i 75 MeVj where we 
set a(0) = 1. 

With the solutions at hand we can backtrack from the 
time of nucleosynthesis (a well-studied period m) to 
the pion gas at temperatures two orders of magnitude 
higher, since the entire particle content in this epoch 
is well known. We do not resort to usual textbook 
power-law approximations since simple computer codes 
produce the (numerically) exact solutions for this one¬ 
dimensional evolution-equation set. For computer accu¬ 
racy, it is necessary to set a unit system that minimizes 
the number of large powers. We often take (100 MeV) 
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Figure 1: Computed scale factor as a function of time (solid, 
blue online). We take as normalization, a(0) — 1, where 
T{t = 0) = 175 MeV. We also show the analytic form 
a{t) oc \/i (dashed line, red online) that fits the lepton era 
at smaller temperatures. As we focuse here on the pion 
gas, found at earlier times, roughly from tr^ns = Os to 
tT^ioo ~ 10“^ peV“^ (10“® s), that square-root approxima¬ 
tion separates significantly from the actual numerical compu¬ 
tation that we employ. 


for temperature, energy and chemical potential {ks = 1) 
and peV“^ for time and space (c = 1). With this, the 
Cavendish constant G turns out to be 1/1.44 (100MeV)^. 
Dimensionally, time is an inverse energy, so that 1 s —>■ 
1.52 X 103peV\ 

The resulting scale factor is shown in figure 

B. Computation of the entropy 

For the calculation of the entropy in the homogeneous 
case we first note that the thermodynamic magnitudes 
are a function of the temperature only. We will also as¬ 
sume thermal equilibrium, vanishing chemical potentials 
and adiabatic expansion. Then the conservation of the 
entropy per co-moving volume implies sa^ = soOq ’"^ti^re 
s = s{T) is the entropy density (see for example [TS] ). 

The second principle of thermodynamics gives the 
total-entropy increase as 

TdS = d{pV) + PdV. ( 8 ) 

where S = sV is the total entropy and V is the volume. 
From this equation it is possible to get the thermody¬ 
namic relations: 

1 , , dP , . 

+ = ( 9 ) 

Therefore the entropy in co-moving volume V is con¬ 
stant and proportional to (t) . Enumerating all the 

species (in equilibrium at the same universal T) 

Pi + ■ ■ ■ + Pn + Pi + ■ ■ ■ + Pn ,, 



Figure 2: (Color online). Number densities of the most 

abundant species during the hadron-lepton epoch. From top 
to bottom: e“, (blue); p~, dashed line (green); 7 
(black); tt (red, solid); v (black, dotted). The number density 
of nucleons is completely negligible during the entire time 
interval. 

The pion gas can be near chemical equilibrium because 
the pion production rate (through 77 -f-)- tttt followed by 
7r“'"7r“ 77 TT^TT^, and similar lepton-lepton inelastic inter¬ 
actions) is sufficient to offset the pion decay rate. In 
figure we show the time-evolution of the number den¬ 
sity for the more relevant temperature span between 175 
and 70 MeV. During this time interval, pions (and also 
muons) are abundant, comparably to the (quasi)massless 
species. 

The entropy density may also be written as oc gs{T) , 
gs(T) being the number of effective degrees of freedom. 
In particular, for ultrarelativistic particles, 

siT)=gs^T\ ( 11 ) 

Due to their relativistic behavior throughout our en¬ 
tire temperature range, the effective number of degrees 
of freedom for photons, electrons and neutrinos is con¬ 
stant. The massive species see drops when T < rrii as 
they become non-relativistic. Our numerical computa¬ 
tion of the entropy can be casted in terms of gsiT) and 
it is plotted in figure In particular, the contribution 
of nucleons (as well as all strange and higher-flavor par¬ 
ticles, not mentioned further) to the entropy density is 
completely negligible. 

III. ENTROPY PRODUCTION BY LOCAL 
DEPARTURES FROM (THERMAL) 
HOMOGENEITY 

A. Solution to the heat equation 

In this section, we consider separations from the homo¬ 
geneous background described in section |T^ For simplic¬ 
ity we will take inhomogeneities to be spherical bulbs 
at temperature different from the background. Thus, 
the temperature field T(r,t) will now also acquire a de¬ 
pendence on position. Local thermal equilibrium as well 


s = 
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Figure 3: Top plot: aggregated effective number of relativistic 
degrees of freedom Qa as a function of /3 = l/T from numer¬ 
ical calculation. Bottom plot: effective number of degrees of 
freedom for pions (solid line, blue online), photons (horizontal 
solid line, black online), electrons and positrons (dotted hori¬ 
zontal line, green online) and muons (dashed line, red online). 
Note that at the highest part of the temperature interval, pi¬ 
ons provide a larger contribution to the entropy density than 
photons, though leptons are the largest carriers of entropy. 


as chemical equilibrium is still assumed (and departures 
thereof can be separately considered in further investiga¬ 
tions that we do not attempt here). 

The departure from the background modifies temper¬ 
ature and entropy density 

T(r,t) = Tback{t) + 5T{v,t), 
s{r,t) = Sback{t) + 6s{r,t). (12) 


Setting as simplest initial condition a bubble of higher 
T than the surroundings, the temperature profile of such 
bulb will evolve according to the heat equation. Then, 


A{ST{r,t)) 


k{T) d{ST{r,t)) 
Cp{T) dt 


(13) 


with k{T) being the heat conductivity. Here the 
constant-pressure specific heat Cp is defined as the deriva¬ 
tive of the background entropy (neglecting the newly pro¬ 



Figure 4: Dotted blue online: Thermal conductivity k as func¬ 
tion of temperature at zero chemical potential from solving 
the Boltzmann equation. Solid (red online): simple interpo¬ 
lating function employed in the heat equation solver. 


duced one) with respect to temperature at constant P: 


Cp{T) 


dSback{T) 


dT 


(14) 


Since we already calculated the contribution of pions to 
the entropy density we can immediately compute 

the partial specific heat of the pion gas (we will further 
drop the superindex tt in this section, as all quantities are 
refered to the pion gas alone). The other non-trivial func¬ 
tion is k(T), the thermal conductivity, which depends on 
the temperature alone and is known from recent and ear¬ 
lier studies. The numeric data [4] from a variational so¬ 
lution of Boltzmann’s equation following the Chapman- 
Enskog expansion is shown in figure Since Cp(T) and 
k{T) are nontrivial functions of the temperature, the heat 
equation does not admit an immediate analytical solu¬ 
tion, so we numerically solve it by brute force with the 
simplest parabolic solver for a partial differential equa¬ 
tion based on the finite-difference method in space and 
the Euler method in time. Thus, in figurej^we also show 
a simple interpolating function for the conductivity in 
the temperature interval of interest that we employ to 
speed up the computer code. 

The valley in the conductivity at mid-temperatures oc¬ 
curs because of the ~ /„■ scales; the dropping low- 
temperature behavior can be obtained from the tttt scat¬ 
tering length and non-relativistic kinetic theory, and at 
high-T dimensional analysis dictates k cx as visible. 
The detailed calculation with the full machinery of phase 
shifts, unitarity, chiral perturbation theory, etc. has been 
reported elsewhere |4]. 

The numeric solution of Eq. ( [T^ , ST{r,t), is shown 
in figure with an initial condition that has a spherical 
profile Gaussian in the radius. 


ST{r,0) = 5Toe~^ . 


(15) 
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Figure 5: Temperature profile T{r,t) of an inhomogeneity of 
initial size R = 2.5 x 10® fm, as a function of the radius r 
for increasing times. Top, solid black line: initial condition 
r(0,0) = 170 MeV. Brown solid line, much flatter of the 
bottom: T{r,tr), tr ~ 10“'^^ s. Other lines illustrate the time 
evolution of the inhomogeneity at intermediate times. 


Here STq is the initial central temperature of the inhomo¬ 
geneity over that of the background, and R is the typical 
radius. 

There are several considerations to choose the size of 
the inhomogeneity. At the largest scale, we can ask our¬ 
selves what is the largest possible radius that will ho¬ 
mogenize during the pion gas lifetime. We must also 
take the size of the bulb small enough so as to respect 
CMB constraints. Glancing back to figure[^ we estimate 
the Hubble horizon reached during the pion gas to be 
around 10“^ — 10“^ peV“^. This means that no homog¬ 
enization can take place over distances larger than about 
a light second (1 — 10) x 10“^ peV“^, or squaring and 
inverting, R must be no larger than « 10^® — lO^^fm. 
This guarantees that the thermal flattening of the bulb 
never violates causality. Further, since the first order 
heat equation is not relativistically causal and we have 
not examined the 2"^^ order formalism, we have to re¬ 
strict ourselves to even significantly smaller spheres. A 
further consideration is that if the inhomogeneity is too 
large, its relaxation time will be so great that when it 
reaches thermal equilibrium, there are no pions left (they 
are abundant for T^ack ~ 175 — 80 MeV). For this reason 
(exclusively of simplicity), we will restrict the study to 
inhomogeneities no bigger than i? « 10® fm. These are 
small enough not to perturb the metric significantly, so 
we can treat them simply as Newtonian perturbations. 
Finally, when we consider the smallest radii of the inho¬ 
mogeneity, in the typical nuclear scale or somewhat more, 
RHIC guidance is available. 


B. Entropy increase in one inhomogeneity 


The variation of the entropy of our inhomogeneity of 
volume V during the relaxation process can be written 
as: 

dSr dSy dSv 

dt dt ^ dt ' ^ ’ 

where St denotes the total entropy, dSy represents the 
entropy exchanged with the rest of the universe and dSy 
the inner entropy production. The exchanged entropy 
dSy can be obtained by means of an integral of the in- 
comming entropy current over the surface of the inhomo¬ 
geneity dV. We will consider the exchange as positive if 
entropy is supplied to the subsystem by the surroundings. 
The entropy current will be denoted by j^- Concern¬ 
ing the internal entropy production dSy we introduce 
the rate of entropy production Gg per unit volume and 
unit time inside the system. In terms of these quantities, 
dSy/dt and dSy /dt may be written as 


dSy 

dt 

dSy 

dt 



■ ndT, , 


tJs dV. 


(17) 


Expressing Eq. (16) in terms of the entropy current 
and density we have: 


dSr 

dt 


d 

dt 


ST dV 


dE js • n-f / dV Gg 


IdV 


(18) 


and use of Gauss’s theorem yields the equation 


dST 

dt 


-V • js + CTs . 


Eor small flows, linear laws hold, such as the Eourier law 
for the heat flux: 


je = -k{T)VT ; (19) 


where je is the heat current vector. Other examples of 
linear laws are Tick’s law for a flavor i concentration flux, 
jj = —DiWrii, with Di being a diffusion coefficient for the 
particle species i; or Ohm’s law for the electric current 
density jg = —k'V4> with jg being the electric current, (/> 
the electric potential and k the electric conductivity. A 
general form for the entropy production Gg is 



Afc Vk 

T 


+ K 


I • Jq 

T 


( 20 ) 

with Ak,Vk the activities and the stoichiometric coeffi¬ 
cients for the ith species involved in inelastic particle re¬ 
actions. In the following we will consider the entropy 
production Gg for the thermal flow alone (first term). 
Basic thermodynamics yields 


dll = T dST = (je ■ n) dE dt 


( 21 ) 
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where U is the internal energy of the inhomogeneity. In¬ 
tegrating over the surface and in time, and using Gauss’s 
theorem, we find the entropy produced in the process of 
relaxation of the inhomogeneity: 


AS't= / / dVdtV- 

Jdv T Jv 


( 22 ) 


Applying Fourier’s law in Eq. (19) we find: 


that is short enough so that the background tempera¬ 
ture does not vary appreciably and can be considered 
constant (T = Tq). That means in particular that k and 
Cp also remain constant (in fact the inhomogeneity has 
not fully spread in this case, but we can deal with this 
numerically later). Then we can make the replacement 


<To) 


as ~ — jST^r'^e 

T2 ^4 0 


2R2 


(29) 


A St 


dVdtV ■ 




Applying now Leibnitz’s rule we get 

ASt = j dVdt (|V5r|2 - TAST) , 


wherein Tq = T^ack + dT^. For a temperature interval 
(23) from 175 MeV to 170 MeV (t G [0 - lO^^] MeV"!), we 
keep k{T)/Tq unchanged and of order one, thus, CTs oc 

-^r^e~ . Carrying out the integral over space, one 

gets 


(24) 


dV a. 


oc 6T^ R. 


(30) 


which is positive, ASt > 0, since AT = A{dT) < 0 
(remember that T^ack was position-independent). 

Comparing with Eq. ( |17[ ) we find the production of 
entropy and the divergence of its flow 


= ■^^|V(5T(r,f)p , (25) 

-V-j,(r,t) = '^AST{r,t). (26) 

The internal entropy produced in dissipating an inho¬ 
mogeneity is an integrated entropy A5'y, obtained from 
the entropy-density production as after integrating over 
the time and space when and where the inhomogeneity 
was relevant, 


ASv{STo) = J dVdtas{r,t). (27) 

To ascertain the size of this produced entropy and assess 
its relative importance, it is natural to quotient it by the 
background entropy in the same volume. Stack, that for 
a spherical disturbance integrating up to the radius R 
(defined above in Eq. (151 as the characteristic Gaussian 
fall-off radius) is 


^back{.dS,Ttack) — -TfR^Sback{Tback) ■ (28) 

We now have all necessary equations and can proceed 
to the numerical computation. 


IV. NUMERICAL RESULTS 
A. One inhomogeneity only 


To put some numbers, take an inhomogeneity of size 
10® fm at STo « 10 MeV; one has then f dVas — 10^ 
MeV. This element multiplied by a time interval At « 
10^® MeV“^ gives an integrated entropy of order 10^®. 
Nevertheless, since at 175 MeV stack is numerically of 
order 10® MeV, the background entropy Stack given by 
Eq. (28) is « 10^^, so the ratio ASy/Stack is « 10“"*. 
Inasmuch as we are considering just a tiny time interval 
in which the bubble did not have enough time to evolve, 
the value of the entropy produced over the entire life of 
the bubble must be larger than this figure, and thus not 
negligible at all (but requires a numeric computation). 

Now, by solving the heat equation for T(r,t) in 
Eq. (|13[ ), we can compute the integral in Eq. (27) with 
Eq. (|25[) and thus numerically obtain ASy- 


Table II: Values of ASy / Stack (in units of 10“®) for different 
temperatures and, for each column, a value of the inhomo¬ 
geneity size given by i?i = 2.5 x 10® , i ?2 = 2.5 x 10^ ,R 3 = 
2.5 X 10® , i ?4 = 2.5 X 10® in fm units. Temperatures are in 
MeV. 


T})ack 

m 

1 

2 

3 

4 


40 

56.6 

46.2 

46.2 

43.6 


35 

46.7 

35.6 

35.6 

33.6 


30 

37.7 

26.3 

26.3 

24.8 


25 

30.1 

18.4 

18.4 

17.4 


20 

24.6 

11.8 

11.8 

11.1 


15 

19.8 

6.7 

6.7 

6.3 


10 

15.9 

3.0 

3.0 

2.8 


5 

18.9 

0.8 

0.8 

0.7 


40 

162.3 

132.7 

132.7 

125.0 


35 

132.6 

101.9 

101.1 

95.3 


30 

105.6 

73.9 

73.9 

70.0 


25 

83.6 

51.0 

51.0 

48.1 


20 

67.5 

32.4 

32.4 

30.6 


15 

53.5 

18.0 

18.0 

17.1 


10 

42.5 

7.9 

7.9 

7.5 


5 

48.7 

2.0 

2.0 

1.9 


To check the computer codes and understand the typ¬ 
ical order of magnitude, let us consider a time period 


Table [IT] shows the numeric computation of ASy 
divided by Stack for different choices of Ttack and STq 
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Figure 6: ASv/Shack against STq for Tiack = 100 MeV and 
different radii R, from top to bottom: R = 2.5 x 10® fm (blue 
online), R = 2.5 x 10® fm (green online), and 2.5 x 10® fm 
(black online). 
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Figure 7: ASv /Shack as a function of R and STo at different 
Thack in the moment of the formation of the inhomogeneity. 
From top to bottom the surfaces correspond to Thack = 100 
MeV (green online) and Thack = 130 MeV (blue online). 


(initial, central intensity of the perturbation). In figure 
we plot the same quantity ASy /Shack against STq 
for different initial sizes. As expected, the bigger the 
bulb is, the more entropy it produces, also relative to 
the background. In figure we simultaneously plot 
ASv/ Shack against the size R and intensity STq of the 
inhomogeneity. It is interesting though expected to note 
that at lower background temperatures the integrated 
entropy becomes larger for equal STq. Mat hem atically 
this comes from the term |VT|® in Eq. (25), which 
increases as the Thack decreases, giving rise to a larger 
entropy production. 


B. Multiple inhomogeneities 

In the early-universe hadronic gas there is no reason 
to think that only one bubble of different temperature 


N-1 N 




2R 4R 2R 4i? 4i? 2R 4i? 2R 

Figure 8: Sketch of the inhomogeneities arrangement. 


would form (as opposed to say a nuclear collision which 
is a system of very limited size). In the absence of data all 
we can give is an upper bound to the entropy produced by 
disposing as many inhomogeneities as possible (as long as 
the background does not lose its meaning). We adopt as 
an extreme limit the density of bubbles when their Gaus¬ 
sian two-sigma walls touch. Thus, we will consider for 
geometric simplicity a Cartesian arrangement featuring 
inhomogeneities disposed as in a simple centered cubic 
structure. The typical size of each inhomogeneity will be 
« 2R in diameter. We take AR as reasonable average sep¬ 
aration between inhomogeneities. The edge of such cube 
has a length of 2R N due to the presence of N inhomo¬ 
geneities plus {N — 1)4 i? due to the {N — 1) spacings, as 
we show in figure]^ 

The background entropy for N inhomogeneities 
occupying a volume Vc is then 

SZi = VcShack, ( 31 ) 

with Va = [(Af - l)4i? -b 2NRf. 

Next we need to model the intensity of each perturba¬ 
tion, JTq. In a plasma this will be randomly distributed. 
Conceivable noise models are white noise (all STq equally 
likely) or Brownian noise (the distribution falls as 1/STg). 
An interesting intermediate case that is ubiquitous in 
physics is the so called 1/f noise [H] that distributes 
the bubbles in proportion to I/STq. Both 1/f and 1 //^ 
noises obviously assign lower density to higher STq We 
currently have no reason to prefer one or another distri¬ 
bution, so we examine all three of them. In future work 
we will examine acoustic oscillations of the gas perform¬ 
ing a spectral analysis so that the coefficient amplitudes 
of each Fourier mode will be left arbitrary, to improve 
the treatment here). 

The noise function is in all three cases of the form 

W) = ^, (32) 

wherein fS = 0,1,2 for the white, 1/f and 1/P noises 
respectively. The normalization constant C is determined 
from the total number of inhomogeneities N in Vc by 



with STa , STh the lower and upper limits respectively 
for the initial temperature of the inhomogeneity, i.e.. 
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Figure 9: against R at fixed Ttack = 100 MeV 

and 5To = 40 MeV. From bottom to top: 1//^ noise function 
(blue online), 1/f noise (red online), constant white noise 
(green online). 


Figure 10: against 6To at R = 2.5 x 10® fm 

and Thack ~ 100 MeV. Blue online: 1//^ noise function, red 
online: 1/f noise, green online: constant white noise. 


6Tq € [STa, 5Tf/\. Too high initial temperatures will in¬ 
volve the quark and gluon plasma and are thus out of 
our reach here, so 5Tb ~ 40 MeV seems reasonable for 
this exploration. As for the smallest dTg taken, since we 
work in the isospin limit (for example in the computa¬ 
tion of the thermal conductivity), it doesn’t make sense 
to retain scales smaller than about 5 MeV where quark- 
mass or electromagnetic isospin breaking effects may play 
a role. Thus, we will choose a temperature interval be¬ 
tween 40 MeV to 5 MeV for the separation above the 
thermal background. In figure we plot 
for different initial sizes at fixed Tback = 100 MeV for the 
three noise profiles. 

The summed, integrated entropy AS^'^ is defined as 
the integrated entropy summed over all inhomogeneities 
weighted by function picked, namely. 


ASl^\5To) = ^ ASv{5To)P{5To) , (34) 


STo 


with ASv defined in Eq. (271. Tables III|IV and 0 give 
the summed, integrated entropy with distribution of in¬ 
homogeneities following the noise functions 1//^ , 1// , 1 
respectively. In the case of the 1/f , 1//^ noises, there is 
a balance between the larger entropy production in the 
hotter bubbles and the larger probability of finding the 
colder ones, yielding a relatively flat entropy-production 
dependence on STq. As it is shown in figure the 
largest entropy production is attained for the white noise 
distribution. Nonetheless, note that in no case AS[^'^ is 
much greater than about 


The outcome of the computation is that the entropy 
produced is only a small fraction of the background en¬ 
tropy in the same volume, because many of the inho¬ 
mogeneities are just of small intensity with these weight 
functions. But these are ad hoc: independent ways of 
assessing what inhomogeneities are possible need to be 
found and we look forward to progress in that respect. 


Table III: (in units of 10 ®) for different tem¬ 

peratures. The four columns of data correspond to i?i = 
2.5 X 10® , R 2 = 2.5 X 10'^, As = 2.5 x 10® , R 4 = 2.5 x 10® in 
fm units. Tback and STq are expressed in MeV. 





Noise function 1//® 


Tback 

STo 

1 

2 

3 

4 


40 

8.1 

5.3 

5.3 

5.0 


35 

7.2 

4.5 

4.5 

4.2 


30 

6.6 

4.1 

4.0 

3.4 


25 

6.3 

2.9 

2.9 

2.6 


20 

4.3 

2.5 

2.5 

2.4 


15 

4.0 

1.5 

1.5 

1.4 


10 

6.2 

9.8 xl0“^ 

9.8 xlO"^ 

9.1 xlO"^ 


5 

7.2 

2.5 xl0“^ 

3.1 xlO"^ 

1.4 xlO"^ 


40 

13.0 

10.6 

10.6 

10.0 


35 

11.9 

9.0 

9.0 

8.5 


30 

10.6 

7.5 

7.5 

7.0 


25 

9.6 

5.8 

5.8 

5.5 


20 

8.9 

4.3 

4.3 

4.0 


15 

8.1 

2.7 

2.7 

2.6 


10 

7.6 

1.4 

1.4 

1.3 


5 

10.5 

4.1 xlO-i 

4.1 xlO"^ 

3.9 xlO-i 


V. CONCLUSIONS 


In this work we have examined the production of en¬ 
tropy by thermal inhomogeneities in the pion gas pro¬ 
duced after the quark-gluon plasma hadronization at the 
early universe. 

In view of the uniformity of the CMB at large scales, 
too large to not have been in causal contact, standard 
theory invokes a time of accelerated expansion (inflation) 
of a universe in thermal equilibrium. Thus, it would ap¬ 
pear natural to assume that such equilibrium was also 
reached early-on at small scales. Nevertheless, the nat- 
urality argument is not fool proof (recall the recent dis¬ 
covery of a light Higgs with no accompanying supersym¬ 
metric partners for any of the SM particles) and does 
not discard the possible existence of small-scale inhomo- 
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Table IV: (in units of 10 ®) for different tem¬ 

peratures. The four columns of data correspond to i?i = 
2.5 X 10® , i?2 = 2.5 X 10^ , R 3 = 2.5 x 10® , R 4 = 2.5 x 10® in 
fm units. Tback and 5Tq are given in MeV. 





Noise function 1// 


Ti)ack 

JTo 

1 

2 

3 

4 


40 

9.8 

7.8 

7.8 

7.0 


35 

7.9 

5.9 

5.9 

4.8 


30 

6.2 

4.4 

4.4 

3.8 

120 

25 

4.9 

3.6 

3.6 

2.2 

20 

4.6 

1.8 

1.8 

1.6 


15 

5.5 

1.1 

1.1 

9.4x10“^ 


10 

4.6 

7.2x10"^ 

7.2x10“^ 

5.3xl0~^ 


5 

5.7 

1.8x10"^ 

1.8x10“^ 

1.4x10“^ 


40 

18.9 

15.4 

15.4 

14.5 


35 

15.6 

11.9 

11.9 

11.2 


30 

12.6 

00 

bo 

8.8 

8.3 

100 

25 

10.0 

6.1 

6.1 

5.8 

20 

8.2 

3.9 

3.9 

3.7 


15 

6.6 

2.2 

2.2 

2.1 


10 

5.3 

9.9x10"^ 

9.9x10“^ 

9.4x10“^ 


5 

6.3 

2.5x10"^ 

2.5x10“^ 

2.3xl0~^ 


Table V: (in units of 10 ®) for different tem¬ 

peratures. The four columns of data correspond to Ri = 
2.5 X 10® , R 2 = 2.5 X 10^ , R 3 = 2.5 x 10® , R 4 = 2.5 x 10® in 
fm units. T^ack and 5Tq are both given in MeV. 





White 

noise 


Tltack 

STo 

1 

2 

3 

4 


40 

195.1 

166.3 

166.3 

149.8 


35 

195.5 

120.1 

120.1 

112.6 


30 

183.4 

120.6 

120.5 

107.1 


25 

153.6 

84.7 

84.7 

78.7 


20 

137.3 

64.3 

64.3 

50.4 


15 

126.1 

30.4 

30.4 

24.3 


10 

57.3 

12.8 

12.82 

12.3 


5 

52.2 

2.7 

2.7 

2.6 


40 

386.5 

315.9 

315.8 

297.5 


35 

315.7 

240.7 

240.7 

226.9 


30 

251.5 

175.9 

175.9 

165.9 


25 

199.0 

121.3 

121.3 

114.5 


20 

160.6 

77.1 

77.0 

72.8 


15 

127.5 

43.0 

42.9 

40.6 


10 

101.1 

18.9 

18.90 

17.9 


5 

118.6 

4.7 

4.6 

4.4 


geneities in the early universe. If any such did not have 


enough time to dissipate before the quark-gluon plasma 
decay cross-over, or was produced during that phase tran¬ 
sition, the entropy production is calculable with modern 
nuclear and particle physics theory. We have exemplified 
with the computation of such entropy increase in a ther¬ 
mal inhomogeneity in the pion gas, but the field is ample 
and much more work is possible. 

We have concentrated on the pion gas at a time inter¬ 
val when pions were main contributors to the universe’s 
entropy, even more than photons, and the ones carry¬ 
ing the largest possible inhomogeneity due to their short 
mean free path. 

Of course, further processes involving non-vanishing 
entropy production, such as flavor or momentum diffu¬ 
sion, are expected to be of equal potential importance. 
We leave them for future work. Here we have remained 
within the realm of small, Newtonian thermal perturba¬ 
tions, and the produced AS'y/S'back is thus a small frac¬ 
tion. In future work we plan to examine the damping 
of acoustic oscillations in this phase by examining sub- 
Jeans modes in the presence of dissipation coefficients, 
and try to put quantitative constraints on the maximum 
size of inhomogeneities that can be dissipated. Looking 
at figure we see that the pion gas can easily reduce a 
thermal perturbation at the 30% level down to the 3% 
level, that is, an order of magnitude. A more detailed 
study beyond this first exploration is granted. Mean¬ 
while, for inhomogeneities of that intensity, production 
of entropy is significant, as can be seen in the various 
tables of subsection IIVBI 

To conclude, we have pointed out that entropy pro¬ 
duction deserves being examined in the hadron-lepton 
phase between 1 and 175 MeV, and we have studied in 
detail one example, that of relaxation of thermal inho¬ 
mogeneities in the pion gas. 
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